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Abstract 

Despite  the  great  amount  of  research  devoted  to  the  visibility  problem,  very  little 
attention  has  been  devoted  so  far  to  measure  or  compare  visibility  data  in  a  way  which 
is  consistent  with  “human-based”  goals  as  of  interest  in  typical  military  applications. 
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measure  distances  in  terms  of  their  visibility  properties.  The  objective  of  this  report  is 
to  introduce  a  functional-theoretic  framework  to  deal  with  the  notion  of  visibility  and 
develop  mathematically  well-defined  and  practically  effective  tools  to  process  visibility 
information,  including  the  comparison,  analysis  and  classification  of  visibility  data.  In 
particular,  a  new  notion  of  difference  of  visibility  distance  is  introduced,  as  a  theoretical 
and  computational  tool  for  the  comparison  of  terrain  models  from  the  point  of  their 
visibility  properties. 
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1  Introduction 


The  visibility  problem  is  a  among  the  most  important  research  topics  in  computational  geom¬ 
etry  and  is  of  interest  in  both  military  applications  such  as  surveillance  and  reconnaissance 
and  civilian  applications  such  as  robotics,  geographic  information  systems  and  computer 
graphics.  Due  to  its  importance,  this  problem  has  been  extensively  investigated  and  several 
efficient  algorithms  have  proposed  to  produce  detailed  terrain  models,  compute  the  visibil¬ 
ity  regions  and  address  viewpoint-based  placement  problems  such  as  the  classical  “museum 
problem”  and  its  generalizations  (see,  for  example,  [1,  13,  17,  18]). 

However,  despite  all  this  effort,  certain  very  relevant  aspects  of  the  notion  of  visibility 
have  not  received  the  proper  amount  of  attention.  In  particular,  there  is  no  satisfactory 
way  to  measure  or  compare  visibility  data  in  a  way  which  is  consistent  with  “human-based” 
goals.  This  is  of  particular  importance  to  effectively  assess,  for  example,  how  well  a  visibility 
model  approximates  the  ground  truth,  or  to  compare  different  visibility  models  and  algo¬ 
rithms.  In  addition,  very  little  effort  has  been  devoted  in  the  scientific  literature  to  examine 
the  functional-analytic  properties  associated  with  the  notion  of  visibility.  Indeed,  in  virtu¬ 
ally  all  studies,  the  visibility  problem  has  been  examined  in  the  context  of  discrete  rather 
than  continuous  data.  While  traditional  methods  from  computer  science  and  computational 
geometry  are  very  effective  to  compute  visibility  structures  and  deal  with  many  practical 
issues  related  to  visibility  computations,  our  point  of  view  is  that,  by  making  use  of  the  con¬ 
tinuous  setting,  one  can  take  advantage  of  the  power  of  functional  analysis  to  gain  a  deeper 
insight  into  the  visibility  problem  and  develop  more  sophisticated  tools  for  the  analysis  and 
processing  of  visibility  data.  . 

The  main  goal  of  this  work  is  to  introduce  a  novel  mathematical  approach  to  deal  with 
the  notion  of  visibility.  By  setting  a  proper  functional-analytic  framework,  we  will  develop 
mathematically  well-defined  and  practically  effective  tools  to  better  process  visibility  infor¬ 
mation,  including  the  comparison,  analysis  and  classification  of  visibility  data  in  a  way  which 
is  consistent  with  “human-based”  goals. 

The  report  will  be  organized  as  follows.  Section  1.1  contains  the  notation  and  definitions 
which  will  be  used  throughout  the  report.  Section  2  sets  the  proper  mathematical  framework 
for  the  description  and  characterization  of  the  viewshed  function  and  sets  the  groundwork  for 
next  section.  Section  3  contains  the  main  results  of  this  work.  Specifically,  a  new  definition 
of  Difference  of  Visibility  Distance  is  introduced  and  validated  on  a  number  of  examples  to 
show  its  consistency  with  human-based  goals.  Using  this  notion,  continuity,  equivalence  and 
approximation  properties  on  a  simple  model  of  univariate  terrain  models  are  examined. 

1.1  Notation  and  Definitions 

We  start  by  introducing  the  basic  definitions  and  notations  which  will  be  used  throughout 
this  report. 
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A  topographic  surface  or  terrain  is  defined  as  the  graph  of  a  bivariate  function  f  :  Q  C 
M2  — >  M,  that  assigns  a  height  or  elevation  f(x )  to  every  point  x  in  a  domain  Q  C  K2.  We 
will  assume  that  Q  =  [0,  l]2  and  that  the  function  /  is  bounded. 

In  practical  situation,  the  values  of  the  function  /  are  known  only  at  a  finite  set  S  C  Q 
of  sample  points.  This  gives  rise  a  discrete  version  of  a  terrain  which  are  known  as  a  Digital 
Elevation  Model  (DEM),  consisting  of  elevation  measures  on  the  discrete  set  of  points  S  C  Q, 
where  S  is  either  a  regular  grid,  or  an  irregularly  scattered  set  of  points.  In  the  first  case, 
we  have  a  Regular  Square  Grid  (RSG),  which  is  obtained  by  partitioning  S  into  equally 
sized  rectangles.  Using  this  model,  the  function  /  is  defined  piecewise  over  each  rectangle, 
e.g.,  using  a  constant  or  bilinear  function.  Each  rectangle  can  be  further  divided  into  two 
triangles,  with  linear  interpolating  functions  defined  on  them.  In  the  second  case,  when  the 
points  in  S  are  irregular  samples,  we  have  a  Triangulate  Irregular  Network  (TIN),  which  is 
defined  by  introducing  a  triangulation  T  of  the  set  S.  Also  in  this  case,  /  is  piecewise  defined, 
for  example  by  using  linear  functions  over  each  triangle  of  T .  Hence  in  both  cases  the  result 
is  a  polyhedral  terrain,  which  is  the  graph  of  a  piecewise  linear  function  and  provides  an 
approximation  of  the  continuous  terrain  function  /.  We  refer  to  [2,  9,  10]  for  additional 
information  about  these  discrete  terrain  models. 


Figure  1:  Examples  of  a  natural  (left)  and  a  urban  (right)  terrain  model  showing  the  visible 
and  not  visible  locations  from  a  viewpoint  denoted  as  the  Observer. 


For  a  fixed  x  G  Q,  we  say  that  the  point  ( y ,  f(y)),  y  G  Q,  belongs  to  the  line  of  sight  of  / 
at  x  if  the  interior  of  the  line  segment  through  ( x ,  f(x))  and  ( y ,  f(y))  lies  strictly  above  the 
surface.  In  this  case,  we  write  that 

(■ y,f(y ))  e  LOS/(x,?/), 

the  point  ( x ,  f(x))  is  called  a  viewpoint  and  y  is  visible  from  x.  It  is  clear  that,  if  ( y ,  f(y))  G 
LOS/(a;,  y),  then  (x,  f(x))  G  LOS f(y,  x).  Hence,  two  points  which  belong  to  the  line  of  sight 
of  each  other  are  mutually  visible. 
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More  generally,  we  can  assign  an  elevation  h  >  0  to  a  viewpoint.  Hence,  for  a  fixed 
elevation  h  >  0,  we  say  that  a  point  ( y ,  f(y ))  belongs  to  the  line  of  sight  of  x  with  elevation 
h,  if  the  interior  of  the  line  segment  through  ( x,f(x )  +  h )  and  ( y,f{y ))  lies  strictly  above 
the  surface.  In  this  case,  we  write: 


(: y,f(y ))  e  LOS f(x,y;h), 


and  the  point  (x,  /(x))  is  a  viewpoint  with  elevation  h. 

For  each  x  G  [0,  l]2,  we  define  the  viewshed  of  f ,  at  the  location  x  with  elevation  h  >  0, 
as  the  set 

V/(x;  h)  =  {y  eQ  :  (y,  f(y))  G  LOS/(x,  y;  h)}  C  Q. 

Hence,  V/(x;  h )  is  the  set- valued  function  consisting  of  the  coordinate  points  of  Q  which  are 
visible  from  the  viewpoint  (x,/(x))  with  elevation  h  (see  Figure  1). 

It  follows  from  this  definition  that  if  h\  <  h-2  then  Vf(x\  hi)  C  Vf(x;  h2). 

Equivalent  Definitions 

For  (xi,z\),  (x2,z2)  G  Q  x  M,  let  L((x\,  z{),  (x2,  z2f)  be  the  line  segment  through  the  points 
(x\,zi)  and  (x2,z2).  We  define  the  viewshed  of  f  at  x  G  Q,  with  elevation  h,  as  the  set 

Vf(x;  h)  =  {y  e  Q  :  Vz  e  xy,  f(z)  <  Lz,  where  (z,  Lz)  G  L((x,  f(x)  +  h),  (■ y ,  f(y)))}. 

The  Line  of  Sight  of  f  with  elevation  h  is 


LOS  f(x,y;h) 


L((x,  f(x)  +  h),  (; y ,  f(y))),  if  y  G  Vf(x;  h); 

oo,  otherwise. 


Other  useful  visibility  structures  for  a  terrain  are  the  horizons ,  describing  the  locations 
on  Q  which  belong  to  the  line  of  sight  at  the  viewpoint  x,  and  block  the  view  of  points  lying 
immediately  beyond  them.  Specifically,  the  local  horizons  of  the  viewpoint  x  are  the  points 
p  G  Q  such  that  (p,  f(p))  G  LOS f(x,p;  h )  and  there  is  no  q  G  Q,  q  ^  p,  such  that  p  G  xq 
and  all  points  in  pq  are  visible  from  x.  The  global  horizons  of  the  viewpoint  x  are  the  points 
p  G  Q  such  that  (p,  f(p))  G  LOS f(x,p;  h )  and  for  every  point  q  G  Q  such  that  p  G  xq,  q  is 
not  visible  from  x.  It  is  clear  that  if  p  is  a  global  horizon  point,  then  it  is  also  a  local  horizon 
point. 

Several  variants  of  the  notion  of  viewshed  are  proposed  in  the  literature  [4,  5,  6,  7,  11]. 
Usually,  this  definition  is  introduced  in  the  discrete  setting  and  it  is  denoted  as  a  function  of 
the  viewpoint  rather  than  the  viewpoint  location  in  Q.  Also  recall  that  there  are  generalized 
or  extended  notions  of  viewsheds  [6,  7],  where  the  visibility  is  restricted  within  a  certain 
distance  from  the  viewpoint,  that  is 

Vf(x,  r;  h)  =  {y  G  Q  :  ( yj(y ))  G  LOS f{x,y;h)  and  | (y,f(y))  ~  (x,f(x)  +  h)\  <  r)}. 
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Another  useful  extension  is  the  situation  where,  for  each  points  y  E  Q,  a  richer  information 
than  a  simple  Boolean  classification  (visible  /  not  visible)  is  provided,  function.  Finally, 
it  is  also  possible  to  let  the  elevation  h  be  a  function  of  the  viewshed  location  x ,  so  that 
h  =  h(x). 

Notice  that,  in  our  notation,  we  make  explicit  the  dependence  on  the  terrain  function  /, 
since  this  functional  relation  will  be  specifically  considered  in  the  following. 


2  Analysis  of  viewshed  function 

Our  first  objective  is  to  examine  the  properties  of  the  viewshed  function  Vf(x;  h )  as  a  function 
of  the  viewpoint  location  x  G  Q  =  [0,  l]2  and  the  elevation  h  >  0.  This  will  set  the 
groundwork  for  the  analysis  of  the  difference  of  visibility  function  which  will  be  introduce  in 
Section  3. 

It  is  clear  that,  for  each  x  G  Q  and  each  h  >  0,  the  set  Vf(x;h )  is  non-empty,  since 
it  always  contains  the  point  (x,f(x)).  It  is  also  easy  to  see  that,  without  some  regularity 
assumptions  on  /,  it  is  virtually  impossible  to  deduce  interesting  or  realistic  statements 
about  the  viewshed  function  Vf.  In  particular,  if  /  is  highly  irregular,  then  Vf{x',h )  might 
consist  of  a  single  point.  For  example,  let  /(f)  =  0  for  rational  f  and  /(f)  =  1  for  irrational 
f.  Then,  for  each  rational  x,  Vf(x;h )  =  {(x,f(x))}  (a  singleton  set),  provided  h  <  1. 
It  is  interesting  to  observe,  however,  that  even  if  /  satisfies  some  regularity  it  does  not 
necessarily  follow  that  V/(x;  h )  is  larger  than  a  singleton.  Consider,  for  example,  the  terrain 

model  fix i,x2)  =  (x?  +  x%)  sin  (  ,  }  )  .  Notice  that  /  is  continuous  at  the  origin,  but 

\  vIda:2  / 

V/((0,  0);  0)  is  the  single  point  {(0,  0,  0)}.  While  both  examples  are  not  very  realistic  models 
of  terrains,  they  are  useful  to  illustrate  that,  for  a  more  realistic  setting,  one  has  to  impose 
some  control  on  the  number  of  discontinuities  of  /  and  some  regularity  on  its  derivative  as 
well. 

Hence,  in  the  following,  we  will  assume  that  /  is  a  piecewise  C 1  function,  where  the 
first  derivative  is  bounded.  These  assumptions  are  consistent  with  most  terrain  models  of 
practical  interest,  including  urban  and  natural  terrains. 

For  simplicity,  in  the  following  we  only  consider  the  situation  where  /  is  a  univariate 
function;  that  is  Q  =  [0, 1].  The  extension  to  the  bivariate  case  is  straightforward  for  many 
of  the  results  reported  below.  The  general  analysis  of  the  bivariate  case  will  be  presented  in 
a  successive  study.  Notice  that,  since  /  is  defined  on  a  compact  domain  Q,  it  follows  from 
our  assumptions  that  f  is  uniformly  bounded  on  Q.  Also,  in  the  following,  we  will  always 
assume  that  the  elevation  h  of  the  viewpoint  satisfies  h  >  0. 
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x  |  x+2d 
x+d 


Figure  2:  When  the  viewpoint  location  moves  from  x  to  x  +  d,  the  function  V(x)  decreases 
abruptly.  This  is  due  to  the  fact  that  the  line  of  sight  has  the  same  slope  as  the  graph  of  f 
on  a  interval  I  of  positive  measure.  As  the  viewpoint  moves  to  x  +  d,  the  whole  interval  I 
which  was  visible  becomes  invisible. 

2.1  Continuity  Properties 

It  is  not  difficult  to  show  that,  even  for  a  regular  terrain  /,  the  set-valued  function  Vf(x;  h ) 
is  very  sensitive  to  changes  in  x  and  h. 

Indeed,  let  f  E  C1  on  Q  and  consider  the  function  V(x)  =  p  (Vf(x;  h)).  To  show  that 
V{x)  is  not  continuous  in  general,  consider  the  example  illustrated  in  Figure  2.  In  this  case, 
because  the  graph  of  /  has  constant  slope  on  an  interval  of  positive  measure,  a  small  change 
of  the  observer  location  x  will  change  abruptly  the  size  of  the  viewshed  Vf(x]  h).  It  is  clear 
that  V (x)  is  even  more  sensitive  with  respect  to  variations  in  x  if  /  is  discontinuous. 

The  same  terrain  /  as  the  one  in  Figure  2  can  be  used  to  show  that  Vf(x]h)  is  not  a 
continuous  function  of  h.  In  fact,  it  suffices  to  fix  the  viewpoint  at  x  +  d  and  vary  h  so  that 
the  interval  /  is  or  is  not  visible. 

However,  as  observed  above,  we  have  that: 


hi  <  h2  =$■  Vf(x]  hi)  C  Vf(x\  h-2). 

Hence  the  function  U{h)  =  p  (Vf(x]  h))  is  monotonic  increasing  and  lower  semicontinuous. 
Indeed,  recall  that  a  function  /  is  lower  semi- continuous  at  x0  if  for  every  e  >  0  there  exists 
a  neighborhood  U  of  xq  such  that  f(x)  >  f(x o)  +  e  for  all  x  E  U.  Equivalently,  this  can  be 
expressed  as 

lirninf  f(x)  >  f{x0). 

X^XQ 


2.2  Set-theoretic  Properties  of  Vf 

In  this  section,  we  make  some  observations  about  the  geometry  of  the  set-valued  function 
Vf(x]  h).  Notice  that  the  assumption  h  >  0  will  play  an  important  role  in  the  following. 
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The  first  observation  is  that  Vf(x:  h )  is  a  set  of  nonzero  measure,  since  there  is  at  least 
one  interval  (containing  the  viewpoint  x)  which  is  contained  in  Vf(x\h).  Indeed,  we  have 
the  following  observation. 

Lemma  2.1.  Let  f  be  piecewise  Cl  on  Q,  with  bounded  derivative  everywhere.  Then  there 
exists  an  interval  containing  x  which  is  contained  in  Vf{x\  h), 

Proof.  Let  s  be  a  regular  point  of  /,  that  is,  /  is  C 1  in  a  neighborhood  of  x,  and  let 
\f\x)\  <  M  for  each  x  E  Q.  It  is  clear  that  x  E  V/(x;  h).  By  the  continuity  of  /,  for  e  >  0 
small  enough,  f(y)  <  f(x)  +  h  for  all  y  E  (x  —  e,x  +  e).  In  addition,  let  e  >  0  be  such  that 
e  <  h/M.  Hence,  for  each  y  E  (x  —  e,x  +  e),  the  line  from  (x,f(x)  +  h )  to  ( y,f{y ))  has 
slope  larger  than  M  (indeed,  h/e  >  M).  This  implies  that,  for  any  such  y  E  (x  —  e,  x  +  e), 

(y,f(y)  e  LOS f(x,y,h). 

If  x  is  a  point  of  discontinuity  of  /',  a  similar  argument  will  hold.  Indeed,  if  a;  is  a 
jump  discontinuity,  then  there  will  be  an  interval  I\  =  [x,x  +  e)  or  I2  —  (x  —  e,  x]  such  that 
(y,  f(y )  e  LOS f(x,y;  h)  for  all  y  E  h  or  y  E  J2.  □ 

Notice  that,  even  for  x  a  regular  point  of  /,  Vf(x;  h)  need  not  be  an  open  or  a  closed  set. 
In  fact,  it  is  easy  to  find  examples  where  Vf(x]  h )  contains  half-open  or  half-closed  intervals. 

The  main  result  of  this  section  states  that,  under  our  assumptions  on  /,  then  Vf(x;h) 
must  be  a  finite  union  of  intervals. 

Proposition  2.2.  a  Let  f  E  Cfl([a,  6]),  where  [a,  b]  is  a  bounded  interval.  Then,  for  each 
x  E  Q  and  h  >  0,  the  set  Vf(x;  h)  is  a  finite  union  of  intervals. 

In  order  to  prove  this  result,  we  first  show  the  following  lemma,  showing  that  the  line 
joining  a  viewpoint  to  the  terrain  at  a  local  horizon  must  be  tangent  to  the  terrain  at  the 
horizon. 

Lemma  2.3.  Let  f  E  C1([a,  b ])  and  suppose  that  p  E  (a,  b )  is  a  local  horizon  of  the  viewpoint 
b  with  elevation  h  >  0.  Then  f(p)  =  ,/(b)t-p^  • 

Proof.  Suppose,  by  contradiction,  that  f'(p)  ^  ■  If  >  f(p)  >  0,  then 

(j>,  f(p)  ^  LOS f(x,p]  h ),  because  the  line  L((x,  f(x)  +  h ),  (p,  f(p)))  will  either  not  intersect 
the  graph  of  /  at  p  or  will  meet  an  obstacle  before  intersecting  the  graph.  The  situation  is 
analogous  if  f(p)  <  <  and  also  in  this  case  ( p,f(p )  ^  LOS f(x,p;h).  On  the 

other  hand,  if  f  (p)  >  >  0,  then  there  is  an  interval  (p,p  +  e)  which  is  all  contained 

in  Vf(b;  h).  In  fact,  there  is  an  e  >  0  such  that,  for  all  p*  E  (p,p  +  e),  the  tangent  line  at 
(p* ,  / (p* )  intersects  the  line  L((b,f(b )  +  h),  (p* ,  / (p* ))).  Hence,  p  is  not  a  local  horizon  of 
b.  Similarly,  if  0  >  f(p)  >  >  o,  then  there  is  an  interval  (p,p  +  e)  which  is  all 

contained  in  Vf(b;  h ),  and  also  in  this  case  p  is  not  a  local  horizon  of  b.  □ 

Notice  that,  if  J  E  [a,  b]  is  an  interval  where  f  E  C1  is  strictly  monotonic  and  y\ ,  y2  G 
J  0  Vf(x]  h ),  it  does  not  necessarily  follow  that  [2/1 , 2/2]  G  h).  In  fact,  there  may  be  a 


point  y  G  (2/1 , 2/2)  where  f'(y )  =  UD  anc[  y  js  a  local  horizon.  Using  this  observation, 

we  deduce  the  following  observation. 

Lemma  2.4.  Let  /  G  C'1([a,6]).  If  2/1, 2/2  G  Jfl  Vf(x;h )  and  f\y)  ^  for  all 

y  G  (2/1, 2/2),  then  [2/1 , 2/2]  e  Vf(x;h). 

Proof.  This  follows  from  Lemma  2.3,  by  observing  that  no  local  horizons  are  possible  in 
the  interval  (2/1,  y2)  •  □  . 

Notice  that  Lemma  2.4  is  satisfied  if,  for  example,  2/1, 2/2  £  J  H  Vf(x;  h )  and  /"(d)  f  0  for 
all  y  G  (dis¬ 
proof  of  Proposition  2.2.  If  /  G  C1([a,6])  and  a;  is  a  viewpoint  with  elevation  h,  one 
can  identify  a  set  S  containing  all  the  local  horizons  of  /  as  follows  as  the  set  of  the  locations 
p  where  f'(p)  =  —  .  If  this  set  is  discrete  then  the  set  Vf(x;  h )  can  be  identified  as 

the  union  of  those  intervals  between  distinct  local  horizons.  On  the  other  hand,  the  set  S 
need  not  be  discrete  since  it  may  contain  an  interval  /  corresponding  to  values  where  f  is 
constant  and  it  satisfies  the  equation  f  (p)  =  —  for  all  p  G  I  C  S.  However,  if  this 

is  the  case,  then  the  interval  /  cannot  be  visible.  Hence,  the  set  Vf(x;  h )  is  again  determined 
by  the  discrete  set  of  points  where  f(p)  =  —  .  □ 

3  Difference  of  Visibility  Distance 

Since  in  both  military  and  civilian  applications  it  is  important  to  compare  different  terrain 
models,  it  is  useful  to  introduce  an  appropriate  notion  of  distance  defined  on  terrains.  Un¬ 
fortunately,  the  usual  mathematical  norms  such  as  the  Lp  or  the  Sobolev  norms,  are  not 
particularly  effective  to  compare  terrains  on  the  basis  of  their  visibility  properties.  In  fact, 
as  further  illustrated  in  the  examples  below,  it  is  easy  to  produce  examples  of  terrains  which 
are  close  with  respect  to  an  Lp  norm,  but  have  completely  different  visibility  properties. 
In  contrast  with  those  norms,  we  will  introduce  a  new  notion  of  distance  which  is  not  de¬ 
fined  directly  on  the  space  of  terrain  models  and  will  be  useful  to  precisely  compare  terrains 
according  to  the  properties  of  the  associate  viewshed  functions.  Specifically,  we  define  the 
Difference  of  Visibility  Distance  between  the  terrains  /  and  g  as: 

BV{f,g;h)  =  [  p,  {Vf(x]  h)  A  Vg{x\  h))  dx,  (1) 

JQ 

where  A  A  B  is  the  symmetric  difference  of  two  sets  (defined  as  {A  fl  Bc )  U  ( Ac  fl  B)),  h  >  0 
is  fixed,  and  p  is  the  Lebesgue  measure.  Notice  that  it  is  easy  to  modify  the  above  definition 
by  using  the  generalized  versions  of  the  Visibility  Region  function  mentioned  in  Section  1.1. 
However  these  generalizations  will  not  be  considered  in  this  report.  Finally,  notice  that  a 
notion  of  Difference  of  Visibility  with  some  similarities  to  the  one  given  above  was  recently 
proposed  in  [15]. 
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For  functions  f,g  satisfying  simple  measurability  conditions,  TV(f,g;  h )  is  a  well-defined 
function  with  range  on  [0, 1]  (since  we  have  chosen  functions  with  support  on  the  normalized 
domain  Q  =  [0,  l]2).  In  particular,  TV(f,g ;  h)  —  0  if  /  and  g  are  associated  with  the  same 
viewshed  function  for  all  values  of  the  viewpoint. 

Notice  that  W(f,g]h)  is  not  properly  a  distance,  in  the  sense  that  it  does  not  satisfy 
the  mathematical  conditions  of  a  metric;  in  particular,  the  triangle  inequality  does  not  hold 
in  general. 

3.1  Examples 

As  a  first  test  to  motivate  and  justify  this  notion  of  distance,  we  will  compute  the  Difference 
of  Visibility  Distance  on  a  number  of  simple  models  of  terrains  to  show  that  this  notion 
will  provide  a  satisfactory  way  of  assessing  the  closeness  of  the  viewsheds  associated  with 
different  terrain  models.  In  the  following,  we  assume  that  the  viewpoint  elevation  is  a  fixed 
constant  0  <  h  <  1. 


(a)  (b)  (c) 


Figure  3:  Examples  of  map  functions:  (a)  Example  1,  (b)  Example  2;  (c)  Example  3. 


Example  1 


This  is  illustrated  in  Figure  3(a). 


Let  /(«„*,)  =  |°  lf  <*'•*>)  €  1°’  1](  and  **’‘*».  This  is  illustrated  in  Figure  3(a). 

V  ;  \l  if  (x1,x2)  e  [0,  l]2  and  X1  =  X2  &  v  ; 

Now,  if  x  =  ( Xi,x2 )  is  chosen  so  that  X\  =  x2,  then  Vf(x;  h )  =  [0,  l]2.  If  X\  <  x2,  then 
Vf(x]h )  =  {(xhX2)  :  xi  <  ^2}  l]2.  Finally,  if  x\  >  x2,  then  Vf(x;  h )  =  {(x-i,x2)  :  x\  > 

x2 }  H[0,  l]2-  Notice  that,  in  the  first  case,  the  Lebesgue  measure  p(Vf(x]  h ))  =  1,  but  in  the 
other  two  cases  p(Vf(x]h ))  =  1/2  (that  is,  unless  the  observer  is  at  the  ridge  x  =  y,  only 
half  of  the  map  is  visible  from  any  location). 
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Let  g(x)  =  0  for  all  x  G  Q.  Notice  that  the  set  Vf(x;  h )  A  Vg(x ;  h )  is  empty  if  X]  =  x2, 
while  it  is  a  half-plane  if  x\  %  x2.  Since  fi(Vf(x]h )  A  Vg(x;h ))  =  1/2  for  a.e.  x  e  Q,  we 
conclude  that: 

TV(/,<?;h)  =  ^ 

This  is  consistent  with  the  intuition  that  g  and  /  are  associated  with  very  different 
viewsheds.  In  fact,  the  value  of  £V(/,  g;h)  =  |  can  be  interpreted  as  the  assessment  that 
using  g  as  an  approximation  model  for  f  we  would  be  wrong  50%  of  the  time. 

Notice  that,  using  the  traditional  Lp  norm,  with  1  <  p  <  oo,  to  measure  the  distance 
between  /  and  g ,  we  obtain  ||/  —  g\\p  =  0  (while  ||/  —  plloo  —  1)-  This  would  not  be  an  useful 
measure  for  the  difference  of  the  visibility  sets  associated  with  /  and  g. 


Example  2 


Let  f(x i,x2)  = 


This  is  illustrated  in  Figure  3(b). 


if  (xi,  x2)  G  [0,  l]2  and  X\  %  x2  +  1/2 
if  (xi,  x2)  G  [0,  l]2  and  X\  =  x2  +  1/2 
Direct  computation  shows  that  /x(V/2(x;  h))  =  7/8  and  /i(V/2(x))  =  1/8,  respectively. 
Again,  letting  g(x)  =0  for  all  x  G  Q  and  proceeding  as  above  we  find  that 


TV(/,  <?;  h) 


17  7  1 

8  8  +  8  8 


7 

32 


Notice  that,  if  the  traditional  Lp  norms  are  used  to  measure  the  distance  between  /  and 
g ,  as  in  Example  1,  we  obtain  ||/  —  g\\p  =  0,  for  1  <  p  <  oo,  and  ||/  —  ^Hoq  =  1. 


Figure  4:  Example  f:  (a)  f\  —  1  in  the  dark  set,  fi  —  0  in  the  white  set,  (b)  f2  —  1  in  the 
dark  set,  f2  =  0in  the  white  set;  (c)  f\  —  f2  on  the  sets  R1  and  R3,  f\  %  f2  on  the  set  R2. 


Example  3 

Let  f{x \,x2)  =  1  —  x\  —  x2,  (xi,x2)  G  Q.  This  is  illustrated  in  Figure  3(c). 
For  each  x  G  Q,  we  have  that  V/(x;  h)  =  [0,  l]2.  Hence  /x(V/(x;  h))  =  1. 
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Letting  g(x)  =  0  for  all  x  G  Q,  we  find  that  Vf(x ;  h)  A  Vg(x ;  h)  is  the  empty  set.  for  each 
x  e  Q.  Hence 

TV(f,g]h)  =  0. 

This  is  the  expected  result  since  /  and  g  are  associated  with  the  same  visibility  regions. 

Notice  that,  if  the  traditional  Lp  norms  are  used  to  measure  the  distance  between  /  and 
g,  we  obtain  \\f  -  g\\p  ±  0.  In  particular,  ||  f  -  g\\i  =  1/6  and  \\f  -  ^Hoo  =  1. 


Example  4 

Let  fi(xl,x2)  = 


1  if  (xiiX2)  G  [0,  l]2  and  X\  <  x2 
0  if  (xi,X2)  G  [0,  l]2  and  X\  >  x2 

1 


Let  f2  be  a  perturbed  version  of  f\ :  f2(xi,x2)  = 


if  (xi,  x2)  G  [0,  l]2  and  Xi  <  x2  —  5 
0  if  (xi,  x2)  G  [0,  l]2  and  X\  >  x2  —  § 
We  refer  to  the  notation  of  Figure  4.  We  observe  that  if  x  G  Rl,  then  Vfx  (x;  h )  A  Vj2  (x;  h ) 
is  empty;  if  x  G  i?2,  then  V/,  (x;  h)  A  l/2(x;  h)  is  i?2;  if  x  G  R3,  then  Vft  (x;  h)  A  V/2(x;  h)  is 
empty.  Hence,  combining  these  observations,  we  have  that 


VV(f1J2)  =  fi(R2)fi(R2)=72. 


where  7  =  g(R2)  is  the  area  of  region  of  perturbation. 


3.2  Extensions 

The  Difference  of  Visibility  T)V(/,  g\  h )  provides  a  global  information  about  the  difference  of 
visibility  properties  of  the  terrains  /  and  g.  A  natural  way  to  introduce  a  “local”  notion  for 
the  difference  of  visibility  distance  is  to  consider  the  quantity: 

'DVv{f,g]h)=  [  n(Vf(x;h)  A  Vg(x]h))  dx,  (2) 

Jv 

where  V  C  Q  is  the  viewpoint  domain.  That  is,  the  viewsheds  of  the  terrains  /  and  g  are 
compared  on  a  subset  of  Q  rather  than  on  the  whole  domain  Q.  More  generally,  it  is  useful 
to  introduce  a  weighted  version  of  the  Difference  of  Visibility  Distance  as 

RWw(f,  g;h)  =  f  g  (V)(x;  h )  A  Vg(x]  h ))  w(x)  dx, 

Jq 

where  w(x)  is  a  positive  weight  function.  It  is  clear  that,  in  the  special  case  where  w(x)  = 
Xv(x),  then  we  recover  the  definition  (2).  In  general,  by  choosing  the  function  w  appro¬ 
priately,  it  is  possible  to  compare  the  viewsheds  of  /  and  g  by  weighting  in  different  ways 
viewpoints  located  at  different  regions. 
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Figure  5:  (a)  Map  function  f ,  (b)  Plot  of  g(Vf(x)) ,  as  a  function  of  the  location  x.  For  this 
example,  a  numerical  computation  gives  Tty(f)  =  0.1475  and,  hence,  the  distance  between  f 
and  the  ground  truth  g  is  TV(f,g)  =  0.8525. 

As  illustrated  in  some  examples  in  Section  3.1,  if  g(x)  is  constant  for  all  x  €  Q,  then 
the  Difference  of  Visibility  Distance  TV(f,  g:  h)  can  be  interpreted  as  a  measure  of  the  close¬ 
ness  of  the  terrain  /  to  a  terrain  which  is  completely  visible.  That  is,  TV(f,g]h)  yields 
the  proportion  of  the  map  which  is  non-visible  from  a  typical  location.  Indeed,  a  simple 
computation  shows  that 


£W,  g\h)  =  1  -  /  //  (Vf(x;  h))  dx 

J[o,i]2 

and,  hence,  TV(/,  g;  h )  =  0  if  and  only  if  the  V/(x;  h)  =  1  for  a.e.  x  G  Q.  This  shows  that 
the  number 

£V(/;  h)  =  TV(f,  f',h)=  [  fi  (Vf(x))  dx 

J[  o,i]2 

describes  a  property  intrinsic  to  the  map  /. 

More  precisely,  this  quantity  yields  the  expected  value  of  the  measure  of  the  visibility 
region  of  f  for  a  randomly  selected  observation  point  (assuming  uniform  distribution  on 
[0,  l]2.  This  is  also  the  probability  that  any  two  randomly  selected  locations  in  the  map  f  are 
in  the  Line  of  Sight  of  each  other  (since  the  size  of  the  map  is  normalized  to  1). 

In  particular  that,  for  the  map  /  from  Example  1,  TV(/;/i)  =  1/2;  for  the  map  /  from 
Example  2,  TV(f;  h)  =  25/32;  for  the  map  /  from  Example  3,  TIV(/;  h )  =  1.  Each  of  these 
numbers  gives  precisely  the  probability  that  a  pair  of  randomly  selected  points  is  the  line  of 
sight  of  each  other.  An  example  of  computation  of  the  number  XV(/;  h)  for  a  more  general 
terrain  model  is  given  in  Figure  5  (the  value  of  TV(/;  h)  for  this  example  is  computed  using 
a  numerical  routine  which  adapts  the  Matlab  functions  available  in  the  Mapping  Toolbox). 
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It  is  clear  that  if  TV(/;  h)  —  1  then  a  single  observer  is  able  to  have  visibility  over  the 
whole  map  /.  More  generally,  if  TV(/;  h )  is  “close”  to  1,  then  most  of  the  map  /  is  visible 
from  a  single  location.  This  suggests  that  TV(f:  h )  could  be  useful  as  a  tool  to  deal  with 
viewpoint-based  placement  problems  concerning  the  distribution  of  observations  points  on  a 
terrain. 


3.3  Analysis  of  the  Difference  of  Visibility  Distance 

In  this  section,  we  examine  the  properties  of  the  Difference  of  Visibility  Distance  TV(f,g]  h ) 
defined  by  (1). 

We  will  only  consider  the  univariate  case,  under  the  assumption  that  the  terrains  f,g 
are  piecewise  polynomial  functions  on  Q  =  [0,1];  that  is,  they  are  polynomials  except  for 
finitely  many  points  ti,  -  -  -  tjv  where  jump  discontinuities  are  allowed.  More  precisely,  for 
m  —  0, 1, . . . ,  we  consider  the  spline  function  spaces 

N 

Sm,N  =  {f  ■  f(t)  =  Xli(t),te  Q}, 

i= 1 

where,  for  each  i,  the  term  p[m\t )  is  a  polynomial  of  degree  m,  and  Notice 

that  we  make  no  regularity  assumptions  at  the  knots  {tl5 . . . ,  t/v}-  In  our  notation  for  Sm^, 
the  index  m  is  associated  with  the  degree  of  the  polynomial  and  the  index  N  to  the  number 
of  knots,  corresponding  to  possible  jump  discontinuities. 

Questions  to  be  investigated  include  the  continuity  of  the  distance  of  visibility  function 
on  these  spaces  and  the  construction  of  “best”  approximations  with  respect  to  the  notion  of 
distance  of  visibility.  For  simplicity,  in  the  following  we  will  consider  mostly  the  space  So, tv; 
that  is,  the  case  of  terrain  models  which  are  piecewise  constant.  This  will  allow  us  to  get  a 
sufficiently  clear  insight  into  the  properties  of  the  Difference  of  Visibility  Distance. 

3.3.1  Continuity. 

The  first  observation  is  about  the  continuity  of  the  Difference  of  Visibility  Distance  on  the 
space  S'o,tv- 

Proposition  3.1.  Let  f,g  E  So,n ■  Given  e  >  0,  there  is  a  S  =  S(e,N )  such  that,  with 
probability  1,  ||/  —  g\ |oo  <  8  implies  that 

W(f,g)<e. 

To  prove  this  result,  we  need  the  following  observation. 
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Figure  6:  A  change  in  f  by  5  on  two  consecutive  intervals  changes  the  size  of  the  visible 

re9ion  bV  tS 


Lemma  3.2.  Let  /,  g  £  S0jn,  i.e.,  f(t)  =  Y^iLi  Q  and  suppose  that  f(x0)  +  h  ^  q, 

/or  all  i  =  1, . . . ,  N.  Then,  given  e  >  0,  there  is  a  S  =  S(e,  N )  such  that  \\f  —  g||oo  <  6  implies 
that 

\p(Vf(x0,  h)  -  n(Vg(x0,  h)  |  <  e. 

Proof.  Since  f(x o)  +  h  ^  q  for  all  i,  the  slope  of  each  line  segment  joining  f(x o)  +  h 
to  f(xi)  satisfies  |  tan 6?^ |  >  A  for  some  A  >  0.  Suppose  that  \\f  —  g\\ oo  <  6(e,N)  =  It 
follows  (see  Figure  6)  that,  for  each  pair  of  consecutive  intervals  and  /l+i ,  the  change  of 
visibility  affects  at  most  a  subinterval  of  size  — (V  =  N2f^e.  <  ff.  This  may  happen  for  at 
most  N/2  locations,  so  that 

N  2e 

HVf(x0,  h )  -  //(V5  (x0,  h)|  <  y  y  =  e-  D 

For  a;0  G  Q  and  h  >  0  being  fixed  values,  let  the  mapping  QXOih  be  defined  by 

®x0,h  '■  f  fJ>(Vf(x0,  h).  (3) 

Lemma  3.2  shows  that,  for  each  x0,  h,  the  mapping  $X0)/l  is  continuous  on  the  space  S0jn, 
except  for  a  subset  of  measure  zero.  Hence  we  can  state  that,  with  probability  1,  the  mapping 
&x0,h  is  continuous  on  the  space  This  also  completes  the  proof  of  Proposition  3.1. 

Similar  computations  can  be  used  to  extend  the  above  results  to  the  space  Sj,7v- 

3.3.2  Equivalence  Properties. 

The  next  group  of  results  are  dealing  with  the  issue  of  establishing  a  notion  of  equivalence 
with  respect  to  the  Difference  of  Visibility  Distance.  In  fact,  as  observed  in  Section  3.1, 
terrains  /  and  g  with  the  same  viewshed  need  not  be  the  same. 
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Definition  3.3.  f  and  g  are  V-equivalent  if,  for  all  x  G  Q  and  for  all  h  >  0,  Vf(x,  h )  = 
Vg{x,h). 

Proposition  3.4.  Let  /,  g  G  S0^.  f  and  g  are  V-equivalent  if  and  only  if  f  =  g  +  c,  where 
c  is  a  constant. 

Proof. 

(4=)  It  is  clear  that  if  the  line  segment  from  f(x)  +  h  to  f(y),  y  G  Q,  meets  no  obstacle 
then  the  same  holds  for  the  line  segment  from  f(x)  +  c  +  h  to  f(y)  +  c,  for  any  c  G  M. 

(=^)  Fix  h  >  0  and  let  Vf(x,h )  =  Vg(x,h).  Suppose  that  f(t )  ^  g(t)  for  some  t  G  Q. 
Hence  we  can  write 

N  N 

f(t)  =  ^2ci{f)xii(t),  g(t)  =  5>M), 

i= 1 

where  q(/)  7^  Ci(g)  for  some  i.  Suppose  that  q+i(^)  >  Ci(g)  and  let  x  =  Xi  +  <5  for  some 
5  >  0.  Notice  that  Vf(x,  h)  fl  /*  =  Vg(x,  h )  D  and,  thus,  it  must  be  that 

Q+i (fl1)  -  Ci{g)  =  Ci+i(f )  -  Ci(f). 

A  similar  argument  holds  if  ci+i(g)  <  Ci(g).  Also,  the  argument  can  be  repeated  for  each 
interval  /*.  This  shows  that,  for  each  i,  ci+i(/)  =  ci+i(g)  +  q(/)  —  Ci(g),  so  that  /  =  g  +  c 
for  some  constant  c.  □ 

Lemma  3.5.  Let  f  G  S0jn  and  suppose  that  Vf(x,  h)  —  1  for  all  x  e  Q.  Then  f  is  constant. 

Proof.  Arguing  by  contradiction,  suppose  that  /  is  not  constant  at  to-  Then  /  has  a 
jump  at  t0-  WLOG,  let  /(tp)  <  /(to).  It  follows  that,  for  x  >  t\ j",  Vf(x,h )  <  1  since  an 
interval  around  t(y  will  not  be  visible  from  x.  This  is  a  contradiction  and,  hence,  /  must  be 
constant.  □ 

If  /  G  Sm>N,  m  >  0,  the  argument  above  does  not  hold  any  more.  However,  we  can  make 
the  following  observation. 

Lemma  3.6.  Let  f  G  Si}n  and  suppose  that  Vf(x,  h)  =  1  for  allx  G  Q.  Then  f  is  continuous 
and  monotonic. 

Proof.  Arguing  by  contradiction,  suppose  that  /  is  not  constant  at  t0.  Then  /  has  a 
jump  at  to-  WLOG,  let  /(tp)  <  /(to)-  It  follows  that,  for  x  >  t/,  Vf(x,h )  <  1  since  an 
interval  around  tf  will  not  be  visible  from  x.  This  is  a  contradiction  and,  hence,  /  must  be 
continuous.  A  similar  argument  shows  that  /  must  be  monotonic.  □ 

We  will  now  show  that  the  mapping  defined  by  (3)  is  injective  on  the  space  S0^jy. 
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Figure  7:  A  shift  of  the  viewpoint  by  5  changes  the  size  of  the  visible  region  by  7 
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Proposition  3.7.  Let  f,gE  So, tv  and  suppose  that 

p(Vf(x,  h ))  =  n(Vg(x,  h ))  for  all  x  G  Q,h  >  0. 

Then  f  and  g  are  V-equivalent.  In  fact,  f  =  g  +  c,  where  c  is  a  constant. 

Proof.  The  first  observation  is  that  p,(Vf(-,h))  has  a  jump  whenever  /  has  a  jump. 
Suppose  that  /  has  a  jump  at  x0  G  Q ,  that  is, 

fix 0)  -  /(To)  =  J- 

WLOG  we  can  assume  that  J  >  0.  For  h  sufficiently  small,  the  sampling  interval  to  the  left 
of  Xq  will  not  be  visible  from  xf ,  but  it  will  be  visible  from  xf.  Also,  if  the  viewpoint  is 
chosen  to  be  xf,  then  all  the  region  Vf(x q  ,  h)  is  contained  in  Vf(x q  ,  h).  Hence  p(Vf(xf ,  h))  — 
n(Vf(x 0  ,  h))  >  jy  (recall  that  the  size  of  a  sampling  interval  is  A). 

Suppose  that  x0  is  a  regular  point.  Since  /  G  So, tv,  this  implies  that  /  is  constant  near 
Xo-  If  the  viewpoint  Xq  is  moved  by  A,  then  all  sightlines  emanating  from  f(x 0)  +  h  will 
be  shifted  accordingly.  In  particular,  for  each  interval  of  visibility  or  non-visibility,  the  size 
of  the  visible  region  is  changed  by  an  amount  linearly  proportional  to  A.  This  shown  that 
p(Vf(x ,  h))  is  a  continuous  function  of  x,  if  a;  is  a  regular  point. 

It  follows  by  the  observations  made  above  that  p{Vf(x,  h)),  as  a  function  of  x,  has  a  jump 
discontinuity  at  x  iff  /  has  a  jump  discontinuity  at  x.  Furthermore,  since  for  each  regular 
point  x, 

lim  fAYfjx  +  h))  -  fj,{Vf  {x,  h))  = 

a^o  A  ’ 

where  c  is  a  constant,  then  p(Vf(x,  h )  is  a  piecewise  linear  function  of  x. 

Finally,  let 

/  Oo)  -  f(x 0)  =  J- 

We  can  consider  a  viewpoint  near  a:(|"  with  elevation  h  sufficiently  small  so  that  “small” 
changes  of  the  viewpoint  location  only  affects  the  visibility  of  the  region  near  Xq.  Then,  if 
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f=0 

0  t\  t-2  t'i  £4  *5  te  —  1 

Figure  8:  Example  of  terrain  f  E  So  o-  There  is  no  approximation  of  f  by  f  E  So  3  such  that 
TV(f,f)<\. 

the  viewpoint  xq  is  shifted  to  the  right  by  5,  it  follows  that  p(Vf(x o+5,  h))—/i(Vf(x 0,  h))  —  J  | 
(see  illustration  in  Figure  7.  This  shows  that,  for  6  and  h  fixed,  the  visibility  change  only 
depends  on  the  jump  height  J. 

Combining  all  these  observations,  it  follows  that  if  f  E  S0jn,  then  n(Vf(x,  h ))  completely 
determine  the  locations  of  the  regular  and  irregular  points  of  /,  as  well  as  its  jump  heights. 
It  follows  that  if  n(Vf(x,  h ))  =  p(Vg(x,  h ))  for  all  x  E  Q,h  >  0,  then  /  and  g  can  only  differ 
by  a  constant.  □ 

Open  Question.  Another  natural  question  is:  under  which  condition  is  the  mapping 
<f>  on  Sq,n  onto?  By  the  proof  of  Proposition  3.7,  we  know  that  if  /  G  So, jv,  then 
is  piecewise  linear  on  Q  x  M=.  Conversely,  suppose  that  &x,h(f)  is  a  piecewise  linear  and 
bounded  function  on  Q.  Does  it  follow  that  there  exists  a  function  /  E  S0,n  such  that 
®x,h{f)  =  /^(V/(x;  h ),  where  x  E  Q,  h  >  0? 

Observation.  Consider  the  mapping  on  So, at  given  by 

$h(f)(x)  =  g(Vf(x]h)), 

where  h  >  0  is  a  fixed  parameter.  Then  there  is  a  constant  c  >  0  such  that 

\\®h{f)\\BV  <  c  ll/Hoo, 

where  BV  denotes  the  Bounded  variation  norm.  In  fact,  the  BV  norm  of  'P/l(/)  is  controlled 
by  the  L°°  norm  of  the  (distributional)  derivative  of  n(Vf(x ;  h))  and  this  is  controlled  by  the 
jump  discontinuities  of  /  G  So, at. 

3.3.3  Approximations  with  respect  to  VV  distance. 

One  main  motivation  for  the  introduction  of  the  Difference  of  Visibility  Distance  is  the  need 
to  compute  approximations  of  terrains  which  are  accurate  with  respect  to  their  visibility 
properties.  Because  TV  is  not  a  norm,  the  approximation  of  a  terrain  in  terms  of  the  TV 
distance  is  quite  different  from  a  projection  into  a  subspace  as  frequently  done  in  the  context 
of  Hilbert  spaces. 

Consider  /  E  S0,n-  We  are  interested  in  computing  approximations  /  of  /,  where  /  E 
S'q,m,  M  <  N,  and  the  approximation  is  understood  in  the  sense  of  approximation  with 
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Figure  9:  A 

respect  to  the  TV  distance.  More  specifically,  our  goal  is  to  develop  an  algorithm  such  that, 
starting  from  /  G  So, at,  a  function  /  G  >So,iv/2  is  computed  such  that 

TV(fJ)=  min  TV(f,g). 

9^0, N/2 

This  would  be  useful,  for  example,  in  the  context  of  terrain  models  with  variable  resolution 
[10,  12,  16],  to  compute  lower  resolution  approximations  of  a  terrain  which  are  optimal  with 
respect  to  their  visibility  properties. 

While  we  are  unable,  at  the  moment,  to  describe  such  an  algorithm,  we  have  the  following 
important  observations. 

The  first  observation  is  that,  given  /  G  S^jv,  in  general  it  is  not  possible  to  find  a 
/  G  S'0,at/2  such  that  TV(f,f )  <  |.  To  show  this  consider  the  example  of  terrain  /  in 
Figure  8  which  takes  only  the  values  0  or  1.  Assume  that  h  • Cl.  If  the  viewpoint  is  located 
on  one  of  the  1-valued  regions,  then  n(Vf(x ;  h ))  =  1/2,  while  if  it  is  located  in  one  of  the 
0-valued,  then  n(Vf(x]  h ))  =  1/6.  If  /  is  approximated  by  /  where  only  the  even  knots  are 
preserved,  then  the  best  one  can  do  is  to  set  /  =  1,  so  that  TV(f,  f )  =  1/2.  No  other  choice 
will  give  a  lower  values  for  TV. 

The  example  shows  that,  in  general,  we  cannot  guarantee  to  produce  an  approximation 
such  that  TV(f,  f )  <  1/2.  In  other  words: 

max  min  TV(f,g)  >  -. 

/SSo.JV  <?£S O.JV/2  2 

The  natural  question  is  whether,  given  any  /  G  Ao.jv,  one  can  find  an  approximation 
/  G  »S'o,tv/2  such  that  TV(f,f )  <  Consider  the  example  of  a  terrain  /  G  *Sb,4  illustrated 
in  Figure  9,  with  a  fixed  viewpoint  x.  It  is  clear  that  taking  the  average  value  of  /  over 
consecutive  subintervals  will  provide  no  control  on  the  Difference  of  Visibility  Distance.  As 
an  alternative  strategy,  one  could  assign  the  highest  value  of  two  consecutive  subintervals, 
as  indicated  in  Figure  9  (bottom).  At  this  point,  in  order  to  ensure  that  visible  regions  are 
still  visible,  one  could  adjust  the  height  of  the  new  partition  intervals,  to  compensate  for 
occlusions  due  to  the  new  sampling. 
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4  Summary  of  Results  and  Future  Work 

The  main  contribution  of  this  report  is  to  introduce  a  novel  notion  of  difference  of  visibility 
distance  which  overcomes  the  limitation  of  traditional  mathematical  norms  and  is  consis¬ 
tent  with  human-based  goals  as  understood  in  application  from  surveillance,  reconnaissance, 
navigation,  asset  emplacement,  robotics,  and  situational  awareness.  This  method  is  espe¬ 
cially  designed  to  compare  terrains  from  the  point  of  view  of  their  visibility  properties  and 
the  examples  described  in  this  report  clearly  demonstrate  the  advantages  of  the  new  Dif¬ 
ference  of  Visibility  Distance  with  respect  to  the  traditional  mathematical  norms  which  are 
frequently  used  to  compare  terrains.  Furthermore,  our  new  approach  opens  up  a  completely 
new  perspective  into  the  study  of  the  visibility  problem,  including  new  notions  of  terrain 
approximation,  representation  and  equivalence  which  are  derived  from  the  new  definition  of 
distance.  In  addition,  our  framework  leads  to  an  intrinsic  characterization  of  terrains  from 
the  point  of  view  of  their  viewshed.  Within  the  limited  scope  of  this  project,  only  simpli¬ 
fied  models  of  terrains  have  been  considered,  mostly  in  the  univariate  case.  In  this  setting, 
we  have  established  conditions  for  terrains  to  be  equivalent  with  respect  to  their  visibility 
properties;  we  have  established  conditions  for  continuity  with  respect  to  the  Difference  of 
Visibility  Distance;  we  have  investigated  a  number  of  properties  dealing  with  the  approx¬ 
imation  of  terrains  with  respect  to  this  distance.  Indeed,  the  analysis  of  the  Difference  of 
Visibility  Distance  presented  in  this  report  provides  only  an  exploratory  investigation  about 
the  properties  of  this  new  notion  of  distance  and  its  applications  to  the  study  of  visibility 
problems.  It  is  clear  that  more  advanced  analytical  tools  will  be  needed  to  deal  with  more 
general  problems  of  approximation  of  terrains  using  this  framework.  Also,  a  new  generation 
of  numerical  algorithms  will  be  needed  to  incorporate  these  ideas  into  effective  methods  for 
visibility  computations. 

Future  work  includes  the  extension  of  a  number  of  one- dimensional  theoretical  results  de¬ 
scribed  to  the  two-dimensional  setting.  This  is  expected  to  be  straightforward  in  most  cases. 
More  importantly,  a  deeper  investigation  will  be  needed  about  the  notion  of  approximation 
with  respect  to  TV  distance,  in  order  to  understand  the  ultimate  potential  of  an  approach 
based  on  this  distance.  In  particular,  it  will  be  very  important  to  examine  the  impact  of  the 
TV  distance  into  existing  strategies  for  visibility  analysis  and  computation  based  on  hierar¬ 
chical  or  multiresolution  terrain  models.  Using  these  methods,  the  terrain  is  modeled  with 
variable  local  resolution,  in  order  to  reduce  the  computational  effort  by  focusing  selectively 
on  regions  of  special  interest.  However,  the  effectiveness  of  this  strategy  depends  heavily 
on  the  ability  to  accurately  approximate  terrains  at  coarser  resolution.  This  is  precisely  the 
situation  where  the  newly  introduced  TV  distance  will  be  able  to  provide  more  accurate  and 
effective  strategies  for  terrain  approximations. 
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